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Abstract. We prove that the D4 root system (equivalently, the set of vertices 
of the regular 24-cell) is not a universally optimal spherical code. We further 
conjecture that there is no universally optimal spherical code of 24 points in S 3 , 
based on numerical computations suggesting that every 5-design consisting of 
24 points in S 3 is in a 3-parameter family (which we describe explicitly, based 
on a construction due to Sali) of deformations of the D4 root system. 



1. Introduction 

In |Cohn and Kumar 07j the authors (building on work by Yudin, Kolushov, and 
Andreev in [Yudin 93, Kolusho v and Yudin 94llKolushov and Yudin 97[|Andreev 961 
lAndreev 97] ) introduce the notion of a universally optimal code in S" l_1 , the unit 
sphere in M™. For a function /: [—1, 1) — > R and a finite set C C S 1 ™ -1 , define the 
energy Ef(C) by 

Ef{C)= £ /«c,c'», 

c.c'ec 

where (c, d) is the usual inner product. We think of / as a potential function, and 
Ef(C) as the potential energy of the configuration C of particles on S 71-1 . Note 
that because each pair of points in C is counted in both orders, Ef(C) is twice the 
potential energy from physics, but of course this constant factor is unimportant. 

A function /: [— 1, 1) — * R is said to be absolutely monotonic if it is smooth and it 
and all its derivatives are nonnegative on [—1, 1). A finite subset Co C is said 

to be universally optimal if E f (C ) < E f (C) for all C C with #C = #C 

and all absolutely monotonic /. We say that Co is an optimal spherical code if 
imax(Co) < imax(C) for all such C, where 

imax(C) := max (c, c') 

cc'ec 

is the cosine of the minimal distance of C. A universally optimal code is automat- 
ically optimal (let /(*) = (! - t)~ N or (1 + t) N for large N). 

In |Cohn and Kumar 07| . linear programming bounds are applied to show that 
many optimal codes are in fact universally optimal. Notably absent from this list 
is the D4 root system Cd 4 , which is expected but not yet proved to be the unique 
optimal code of size 24 in S 3 . This root system can also be described as the vertices 
of the regular 24-cell. 
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It is shown in [Cohn and Kumar 07] that the vertices of any regular polytope 
whose faces are simplices form a universally optimal spherical code. The dodeca- 
hedron, 120-cell, and cubes in R" with n > 3 are not even optimal spherical codes 
(see [Sloane 00] ) and hence cannot be universally optimal. Thus the 24-cell was 
the only remaining regular polytope. Cohn and Kumar conjectured in an early 
draft of [Cohn and Kumar 07] that Cd 4 was universally optimal, but reported the 
numerical result that for the natural potential function f(t) = (1 — t)~ x there was 
another code C C S 3 with #C = 24 at which Ef has a local minimum only slightly 
larger than Ef(Coi) (668.1902+, compared with 668). 

What makes this code noteworthy is that in simulations of particle dynamics on 
S 3 under the potential function / (along with a viscosity force to remove kinetic 
energy and cause convergence to a local minimum for Ef), 24 particles converge 
more than 90% of the time to C, rather than to Cx> 4 . Similar effects appear to 
occur for fit) = (1 — t)~ s for other values of s. In other words, these codes have a 
much larger basin of attraction than Cx> 4 , despite being suboptimal. 

In this paper we give a simple description of a one-parameter family of configura- 
tions Cg that includes these codes, and exhibit choices of / (such as fit) = (1 + <) 8 ) 
and 9 for which Ef(Cg) < Ef(CD 4 )- We thus disprove the conjectured universal 
optimality of Cd 4 ■ 

We further conjecture that there is no universally optimal spherical code of 24 
points in S 3 . Any such code would have to be a 5-design, because Cd 4 is. Numerical 
computations led us to a 3-parameter family of such designs that can be constructed 
using an approach introduced by Sali in [Sali 94] . The family contains Cd 4 as a 
special case, and consists of deformations of Cd 4 - 

We exhibit these designs and prove that, within the family, Cd 4 minimizes the 
energy for every absolutely monotonic potential function, and is the unique mini- 
mizer unless that function is a polynomial of degree at most 5. Our computations 
suggest that every 5-design of 24 points in S 3 is in the new family. If true, this 
would imply the nonexistence of a universally optimal design of this size in S 3 
because we already know that Cd 4 is not universally optimal. 

One way to think about the D4 root system's lack of universal optimality is 
that it explains how D4 is worse than Eg- The D4 and Eg root systems are 
similar in many ways: they are both beautiful, highly symmetrical configurations 
that seem to be the unique optimal spherical codes of their sizes and dimensions. 
However, one striking difference is that linear programming bounds prove this 
optimality and uniqueness for Eg but not for D4 (see [Arestov an d Babc nko 971 
IBannai and Sloane 81[|Levenshtein 791|Odlyzko and Sloane 79 ). This leads one to 
wonder what causes that difference. Is D4 in some way worse than Eg? Our results 
in this paper show that the answer is yes: for Eg, linear programming bounds prove 
universal optimality (see |Cohn a nd Ku mar 07] ). while for D4 universal optimality 
is not merely unproved but in fact false. 



2. The codes Cg 

We computed the 24 x 24 Gram matrix of inner products between the points of 
the suboptimal but locally optimal configuration mentioned above for the potential 
function fit) = (1-t)" 1 . Each inner product occurred more than once, suggesting 
that the configuration had some symmetry. By studying this pattern we eventually 
identified the configuration with a code in the following family of 24-point codes 
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Ce C S 3 . (We are of course not the first to use this approach of computing a 
code numerically and using its Gram matrix to detect symmetries and then find 
good coordinates. One recent case — also, as it happens, for codes in S 3 — is 
Sloa ne et al. 03] . where the method is called "beautification.") 
We identify M 4 with the complex vector space C 2 so that 

S 3 = {(w u w 2 ) e C 2 : K| 2 + \w 2 \ 2 = 1}. 

For 9 E R/27rZ such that sin 29 ^ and sin 6* 7^ cos# we set 

Ce :— {(z,0), (0,io), (zsin9,w cos9), (z cos#, w sin 6*) : z 3 = w 3 = 1}. 

Thus Ce consists of 24 unit vectors, namely 3 + 3 of the form (z, 0) or (0,w) and 
3 2 + 3 2 of the form (z sin 9, w cos 9) or (z cos 9, w sin 9) . Each of these codes has 72 
symmetries (each complex coordinate can be independently conjugated or multi- 
plied by cube roots of unity, and the two coordinates may be switched), forming 
a group G isomorphic to the wreath product of the symmetric group S3 with S 2 - 
This group does not act transitively: there are two orbits, one consisting of the six 
points (z,0) and (0,u>) and the other consisting of the remaining 18 points. 

Listing all possible pairs c, d S Cg with c 7^ c', we find that there are in general 
11 possible inner products, with multiplicities ranging from 18 to 84. We thus 
compute that 

E f (Cg) = 18(/(O) + /(sin20)) + 36(/(sin 0) + /(cos 0)) 

+ 36 ^/ ^sin 2 - i cos 2 9^j + f (cos 2 6 - ^ sin 2 I 

+ 72(/(-2f) + /(-^) + /(4^) + /(-^ 
+ 84/(-i 
For (7d 4 we have the simpler formula 

E f {C Dl ) = 24/(-l) + 192 (f + / + 144/(0). 

3. Failure of Universal Optimality 

By Theorem 9b in |Widder "ill p. 154], an absolutely monotonic function on 
[—1,1) can be approximated, uniformly on compact subsets, by nonnegative lin- 
ear combinations of the absolutely monotonic functions f(t) = (1 + t) k with k E 
{0,1,2,...}. To test universal optimality of some spherical code Cq it is thus 
enough to test whether E f {C ) < E f (C) holds for all C C S" 1 " 1 with #C = #C* 
and each f(t) = (1 +t) k . We wrote a computer program to compute Ef(CD 4 ) and 
Ef(Cg), and plotted the difference Ej(Cd 4 ) ~ Ef(Cg) as a function of 9. 

For k < 2 the plots suggested that S/(Cd 4 ) = Ef(Cg) for all 9. This is easy 
to prove, either directly from the formulas or more nicely by observing that C_d 4 
and Ce are both spherical 2-designs (the latter because G acts irreducibly on R 4 ), 
so E f (C Di ) + 24/(1) and E f {C e ) + 24/(1) both equal 24 times the average of 
c 1—* /((c, Co)) over S 3 for any cq E S 3 . 

For k = 3, the plot suggested that S/(Cd 4 ) < Ef(Cg), with equality at a unique 
value of 9 in [0, it], numerically 9 = 2.51674+. We verified this by using the rational 
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paramctrization 

• a 2u r 1 

sine/ = ?r, cost 



1 + V? ' 1 + u 2 

of the unit circle, computing Ef(Cjj 4 ) — Ef(Cg) symbolically as a rational function 
of u, and factoring this function. We found that 

F (C \ F (C \ 1 o (" 6 -6^ 4 -12^ + 3 M 2 -2) 2 
E f {C Di ) - E f (C e ) = -18 — 2 - — 6 , 

and thus that Ef(Cn 4 ) < Ef(Cg), with equality if and only if u is a root of the 
sextic u 6 — 6?i 4 — 12m 3 + 3u 2 — 2. This sextic has two real roots, 

u = -(0.51171+), u = 3.09594-, 

which yield the two permutations of {sin6>, cos9} = {0. 58498-, -(0.81105+)} and 
thus give rise to a unique code Cg with Ef(Co 4 ) — Ef(Cg). This code is char- 
acterized more simply by the condition that sin 9 + cos 9 is a root of the cubic 
3y 3 — 9y — 2 = 0, or better yet that sin 3 9 + cos 3 9 = — |. The latter formulation 
also lets us show that this is the unique Cg that is a spherical 3-design: the cubics 
on R 4 invariant under G are the multiples of Re(w 3 ) + Rc(wf ), and the sum of this 
cubic over Cg is 6 + 18 (sin 3 9 + cos 3 9). Since Cd 4 and this particular Cg are both 
3-designs, they automatically minimize the energy for any potential function that 
is a polynomial of degree at most 3. We must thus try k > 3 if we are to show that 
Cd 4 is not universally optimal. 

For k = 4 through k = 7 the plot indicated that Ef(Cg) comes near Ef(CD 4 ) for 
9 ;=y 2.52 but stays safely above Ef(CD 4 ) for all 9, which is easily proved using the 
rational parametrization. (For k = 4 and k — 5 we could also have seen that Cd 4 
minimizes Ef by noting that Cd 4 is a 5-dcsign.) However, for k — 8 the minimum 
value of E f (Ce), occurring at 6 = 2.529367746+, is 5064.9533+, slightly but clearly 
smaller than Ef(Co 4 ) = 5065.5. That is, this Cg is a better code than Cd 4 for the 
potential function (1 + t) 8 , so Cd 4 is not optimal for this potential function and 
hence not universally optimal. 

The maximum value of E / (C £> 4 ) — E '/ (Cg) for f(t) = (l + t) k remains positive for 
k = 9, 10, 11, 12, 13, attained at values of 9 that slowly increase from 9 = 2.52937— 
for k = 8 to 9 = 2.54122— for k = 13. Each of these is itself enough to disprove 
the conjecture that Cd 4 is universally optimal. (Another natural counterexample 
is f{t) = e 6t with 9 = 2.53719+.) 

We found no further solutions of Ef(Cr> 4 ) > Ef(Cg) with k > 13. It is clear that 
Ef(Co 4 ) < Ef(Cg) must hold for all 9 if k is large enough, because t max (Cg) > 
imax(Cz5 4 ) = \- the smallest value t of t max (Cg) is (V7 - l)/3 = 0. 54858+, 
occurring when either 



with 9 = 2.56092+ or 



to = sin 9 = cos 2 9 — — sin 2 ( 
to = cos 9 = sin 2 9 — ^ cos 2 I 



with 9 = 5.29305+. Quantifying what "large enough" means, and combining the 
resulting bound with our computations for smaller k, we obtained the following 
result: 

Proposition 3.1. For 8 < k < 13, there exists a choice of 9 for which 

E f (Cg) < E f (C D4 ) 
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when f(t) = (1 + t) k . For other nonnegative integers k, no such 8 exists. 
Proof. If k is large enough that 

i8/ (^~^) > 24 /(- i )+ 192 (/(£) +i44/ (°)- 

then Ef(CD 4 ) < Ef(Ce). This criterion is wasteful, but we have no need of sharper 
inequalities. Calculation shows that this criterion holds for all k > 75. This leaves 
only finitely many values of A:. For each of them, we use the rational parametrization 
of the unit circle to translate the statement of Proposition 13.11 into the assertion of 
existence or nonexistence of a real solution of a polynomial in Z[u], which can be 
confirmed algorithmically using Sturm's theorem. Doing so in each case completes 
the proof. □ 

We cannot rule out the possibility that there exists a universally optimal 24-point 
code in S* 3 , but it seems exceedingly unlikely. If C'rj 4 is the unique optimal spherical 
code, as is widely believed, then no universally optimal code can exist. The same 
conclusion follows from the conjecture in the next section. 

It is still natural to ask which configuration minimizes each absolutely monotonic 
potential function. We are unaware of any case in which another code beats Cd 4 
and all the codes Ce for some absolutely monotonic potential function, but given 
the subtlety of this area we are not in a position to make conjectures confidently. 

4. New Spherical 5-Designs 

Spherical designs are an important source of minimal-energy configurations: a 
spherical r-design automatically minimizes the potential energy for f(t) = (1 + t) k 
with k < t. Conversely, if an TV-point spherical r-design exists in S n , then every 
iV-point configuration in S 71 ^ 1 that minimizes the potential function f(t) — (l + t) T 
must be a r-design. Thus, when searching for universally optimal configurations, 
it is important to study r-designs with r as large as possible. 

For 24 points in S 3 , the D4 root system forms a 5-design. By Theorem 5.11 in 
[Delsarte et al. 77] . every 6-design must have at least 30 points, so 24 points cannot 
form a 6-design. Counting degrees of freedom suggests that 24-point 4-designs are 
plentiful, but 5-designs exist only for subtler reasons. One can search for them by 
having a computer minimize potential energy for f(t) = (1 + 1) 5 . Here, we report 
on a three-dimensional family of 5-designs found by this method. The D4 root 
system is contained in this family, and all the designs in the family can be viewed 
as deformations of Cd 4 - We conjecture that there are no other 24-point spherical 
5-designs in S* 3 . We shall show that this conjecture implies the nonexistence of a 
universally optimal 24-point code in S 3 . 

Our construction of 5-designs slightly generalizes a construction of Sloane, Hardin, 
and Cara for the 24-cell (Construction 1 and Theorem 1 in [Sloane et al. 03] V The 
Sloanc-Hardin-Cara construction also works for certain other dimensions and num- 
bers of points, and can be further generalized using our more abstract approach. 
For example, one can construct a family of designs in S* 2 ™ -1 from a design in CP" -1 . 
We plan to treat further applications in a future paper. See also the final paragraph 
of this section. 

Fix an "Eisenstein structure" on the D4 root lattice, that is, an action of Z[r], 
where r = e 27 ™/ 3 is a cube root of unity. (It is enough to specify the action of r, 
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which can be any element of order 3 in Aut(£>4) that acts on R with no nonzero 
fixed points; such elements constitute a single conjugacy class in Aut(£>4).) Then 
R is identified with C 2 , with the inner product given by 

((21,^2), (Ci»C2)) = Re(ziCi + 22(2)- 

The group of sixth roots of unity (generated by —1 and r) acts on Co A and parti- 
tions its 24 points into four hexagons centered at the origin; call them H , Hi,H 2 , H 3 . 
In coordinates, we may take 

H = {(w,0) : w S fi 6 }, 
Hi = {(uiw, tiw) : w £ fie}, 
H 2 = {(uiw, rtiw) : w £ /j,q}, 
H 3 = {(uiw, rtiw) : w £ fig,}, 

where u = ^/l/3 and t — For any complex numbers ao, ai, a%, a 3 with 

|ao| = I ai I = |oa| = \a 3 \ = 1, we define 

D(a , a\,a 2 ,a 3 ) = a H U aiHi U a 2 H 2 U a 3 H 3 . 

We claim that D(a , a\,a 2 ,a 3 ) is a 5-design. That is, we claim that for every poly- 
nomial P of degree at most 5 on R 4 , the average of P(c) over c £ D(a , a\,a 2 ,a 3 ) 
equals the average of P(c) over c £ S 3 . It is sufficient to prove that the average 
is independent of the choice of ao, a-i, a 2 , a 3 , because D(l, 1, 1, 1) = Cd 4 is already 
known to be a 5-design. But this is easy: for each m £ {0,1,2,3}, the restriction 
of P to the plane spanned by H m is again a polynomial of degree at most 5, and 
a-rnHm is a. 5-design in the unit circle of this plane, so the average of P over a m H m 
is the average of P over this unit circle, independent of the choice of a m . 

This construction via rotating hexagons is a special case of Lemma 2.3 in [Sali 94 , 
where that idea is applied to prove that many spherical designs are not rigid. Sali 
rotates a single hexagon to prove that Cd 4 is not rigid, but he does not attempt a 
complete classification of the 24-point 5-designs. 

It is far from obvious that there is no other way to perturb the D4 root system 
to form a 5-design. For example, if there were two disjoint hexagons in D4 that 
did not come from the same choice of Eisenstein structure as above, then rotating 
them independently would produce 5-designs not in our family. However, one can 
check via a counting argument that every pair of disjoint hexagons does indeed 
come from some common Eisenstein structure. This supports our conjecture that 
there are no other 24-point spherical 5-designs in S 3 . 

The family of 5-designs of the form D(ao, 01, a 2 , a 3 ) is three-dimensional, for the 
following reason. One of the four parameters ao,ai,a2,a3 is redundant, because 
for every a £ C* with \a\ — 1 we have D(aao, crai, aa 2 , aa 3 ) = D(ao, a\, a 2 , a 3 ). 
We may thus assume a = 1- We claim that for each (ai,a 2 ,a 3 ) there are only 
finitely many (a'^a^a^) such that D(l, ai, a 2 , a 3 ) = D(l, a^, a^, 03). If this were 
not true, there would be an infinite set of designs D(l, b±, b 2 , b 3 ) equivalent under 
automorphisms of S 3 that stabilize Ho pointwise. But this is impossible, because 
such an automorphism must act trivially on the first coordinate z\. Hence our 
5-designs constitute a three-dimensional family, as claimed. 

Some other known spherical designs can be similarly generalized. For instance, 
the 7-design of 48 points in S 3 , obtained in [Sloane et al. 03] from two copies of Cd 4 , 
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has a decomposition into six regular octagons, which can be rotated independently 
to yield a five-dimensional family of 7-designs. 

It is also fruitful to take a more abstract approach. A 24-point design with 
a fie action is specified by four points, one in each orbit. Our new designs are 
characterized by the condition that under the natural map C 2 \ {(0,0)} — > CP 1 
given by (^1,^2) Z1/Z2, the four points must map to the vertices of a regular 
tetrahedron (if we identify CP 1 with S 2 via stereographic projection, with CP 1 = 
CU{oo} and S 2 a unit sphere centered at the origin). In slightly different language, 
we have specified the image of the design under the Hopf map S 3 — > S 2 . The fact 
that the regular tetrahedron is a spherical 2-design in S 2 plays a crucial role, and can 
be used to prove that this construction yields 5-designs. Likewise, the six octagons 
that make up each of our 7-designs map to the vertices of a regular octahedron, 
which is a spherical 3-design. Again, we intend to discuss this approach in more 
detail in a future paper. 



In this section we prove that among all these new 5-designs, the 24-cell minimizes 
potential energy for each absolutely monotonic potential function. As before, it is 
sufficient to do this for f(t) = (1 + t) . For k < 5 this follows immediately from 
the spherical design property, and for k > 5 we will show directly that Cd 4 is the 
unique minimizer. 

Within each hexagon H m , the six points are in the same relative position in each 
design, and thus make the same contribution to the potential energy. Hence it suf- 
fices to show that the potential energy between each pair of hexagons is separately 
minimized for the D4 configuration. 

Let do = 1, cii = ie 10 , a 2 = ie 1 ^ , and 0,3 — ie 1 ^ . Because of the sixfold rotational 
symmetry of each H m , the angles 9, (f>, and ip are determined only modulo 7r/3. In 
particular, 9 = <fi = tp = tt/6 yields the 24-cell (because tt/2 = tt/6 (mod 7r/3)). 

First, consider the pair (Ho, -Hi)- We find that the inner products between the 
points of Hq and a\H\ are (l/\/3) cos(# + j7r/3) with < j < 5, each repeated 
six times. By Lemma 15.11 below, the sum is minimized exactly when 9 = tt/6 
(mod 7r/3). Similarly, considering (Ho,H2) and (Ho,!!^) shows that the corre- 
sponding contributions to potential energy are minimized when = 7r/6 (mod 7r/3) 
or ip = tt/6 (mod tt/3), respectively. 

Next, consider the pair (Hi, -#2)- The dot products possible are of the form 



with < j < 5, each repeated six times. Once again we conclude from Lemma 
15.11 that the potential energy is minimized when 9 = (j) (mod 7r/3). Similarly, 
considering the remaining pairs (Hi,H^) and {H2,H^) shows that 9 = <p = ip 
(mod tt/3). Thus, it will follow from the lemma below that for each k > 6, the 
D4 configuration with # = = ?/> = 7r/6is the unique code in this family that 
minimizes the potential energy under the potential function f(t) = (1 + t) k . 

Lemma 5.1. Let k be a nonnegative integer. When k > 6, the function 
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has a unique global minimum within [0, ir/3], which occurs at 9 = it/6. When k < 5, 
the function is constant. 

Proof. We must show that the coefficient of y k in the generating function 




j=0 k=0 

is zero if k < 5 or 9 = ir/6 (mod w/3) and strictly positive otherwise. Explicit 
computation using the sum of a geometric series shows that the generating function 
equals 

y 6 (cos 2 (9) (4 cos 2 (9 -3) 2 ( 1 2 2 \ A 1 



216 \l-y "2-y 2-3y 



The factor of (cos 2 9) (4 cos 2 9 — 3) vanishes iff 6 = w/6 (mod tt/3) and is positive 
otherwise. Clearly the factor 

12 2 
T^~7j + Y~y + 2 - 3y 
has positive coefficients, as does 

5 

n 



i 



because l + cos(0+j7r/3)/\/3 > for all j. It follows that their product has positive 
coefficients, and taking the factor of y e into account completes the proof. □ 

6. Local Optimality 

So far, we have not addressed the question of whether our new codes are actually 
local minima for energy. Of course that is not needed for our main result, because 
they improve on the 24-cell regardless of whether they are locally optimal, but it is 
an interesting question in its own right. 

For the codes Cg this question appears subtle, and we do not resolve it com- 
pletely. To see the issues involved, consider the case of f(t) = (1 — As 9 
varies, the lowest energy obtained is 668.1920+ when 9 = 2.5371+. That code 
appears to be locally minimal among all codes, based on diagonalizing the Hessian 
matrix numerically, but we have not proved it. By contrast, the other two local 
minima within the family Cg (with energy 721.7796+ at 9 — —(2.0231+) and en- 
ergy 926.3218+ at 9 — 0.5320+) are critical points but definitely not local minima 
among all codes; the Hessians have 22 and 36 negative eigenvalues, respectively. 

We do not know a simple criterion that predicts whether a local minimum among 
the codes Cg as 9 varies will prove to be a local minimum among all codes, but 
it is not hard to prove that every critical point in the restricted setting is also an 
unrestricted critical point. Specifically, a short calculation shows that for every 
code Cg and every smooth potential function, the gradient of potential energy on 
the space of all configurations lies in the tangent space of the subspace consisting 
of all the codes Cg. It follows immediately that if the derivative with respect to 
6 of potential energy vanishes, then the gradient vanishes as well. Furthermore, 
such critical points always exist: starting at an arbitrary code Cg and performing 
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gradient descent will never leave the space of such codes and will always end at a 
critical point. 

At this point one may wonder whether it is even clear that the regular 24- 
cell is a local minimum for all absolutely monotonic potential functions. It is 
straightforward to show that it is a critical point, but we know of no simple proof 
that it is actually a local minimum. The best proof we have found is the following 
calculation. 

For each of the 24 points, choose an orthonormal basis of the tangent space to 
S 3 at that point, and compute the Hessian matrix of potential energy with respect 
to these coordinates. Its eigenvalues depend on the potential function, but the 
corresponding eigenspaccs do not. There is a simple reason for that, although we 
will not require this machinery. Consider the space Sym 24 (5' 3 ) of all unordered sets 
of 24 points in S 3 . The symmetry group of the 24-cell acts on the tangent space to 
Sym 24 (S' 3 ) at the point corresponding to the 24-cell, and this representation breaks 
up as a direct sum of irreducible representations. On each nontrivial irreducible 
representation the Hessian has a single eigenvalue, and these subspaces do not 
depend on the potential function. In practice, the simplest way to calculate the 
eigenspaces is not to use representation theory, but rather to find them for one 
potential function and then verify that they are always eigenspaces. 

If the potential function is /: [— 1, 1) — > K, then the eigenvalues of the Hessian 
are 



0, 

V"(\) +8/"(0) + 2/"(-i) -12/' Q) +12/' (-^ 



2/"Q) +4/"(0) + 6/"(-£) -8/'Q) -4/'(0) + 8/'(-£) +4/'(-l), 
5f "{l) +4/"(0) + 3/"(-l) -14/'Q) +8 /'(O) + 2/'(-0 +4/'(-l), 

6rQ) + 6.r(-i)-i 2/ '(l) + i2/(-l), 

+ 4 /"(°) + 6 /"(4) +4/ '(^) + 8 /'(°) + 2 °/'(-0 +4/'(-l), 
6/"(i) + 8.r(0) + 6r(-i)-4.f(i) +4 /'(-i), 
8f "{l) + 4 /"( O )- 8 /'(V2)-4/'(O)+8/'(-0 +4/'(-l), 

with multiplicities 6, 9, 16, 8, 12, 4, 9, and 8, respectively. 

One mild subtlety is that is always an eigenvalue, so one might worry that the 
second derivative test is inconclusive. However, note that the potential energy is 
invariant under the action of the 6-dimensional Lie group 0(4), which yields the 
6 eigenvalues of 0. In such a case, if all other eigenvalues are positive, then local 
minimality still holds, for the following reason. Notice that 0(4) acts freely on 
the space of ordered 24-tuples of points in S 3 that span R 4 , and it acts properly 
since 0(4) is compact. The quotient space is therefore a smooth manifold, and the 
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positivity of the remaining eigenvalues suffices for the potential energy to have a 
strict local minimum on the quotient space. 

To complete the proof, we need only consider f(t) = (l+t) k with k G {0, 1,2,...}. 
For k < 5 the other eigenvalues are not all positive (some vanish), but because the 
24-cell is a spherical 5-design it is automatically a global minimum for these ener- 
gies. For k > 6 one can check that all the other eigenvalues are positive. That is 
obvious asymptotically, because they grow exponentially as functions of A;; to prove 
it for all k > 6 one reduces the problem to a finite number of cases and checks each 
of them. It follows that the regular 24-cell locally minimizes potential energy for 
each absolutely monotonic potential function, and it is furthermore a strict local 
minimum (modulo orthogonal transformations) unless the potential function is a 
polynomial of degree at most 5. 
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